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ABSTRACT 

An extremely cold and big spot in the WMAP 1-year data is analyzed. Our work is 
a continuation of a previous paper (Vielva et al. 2004) where non-Gaussianity was 
detected, with a method based on the Spherical Mexican Hat Wavelet (SMHW) tech- 
nique. We study the spots at different thresholds on the SMHW coefficient maps, 
considering six estimators, namely number of maxima, number of minima, number of 
hot and cold spots, and number of pixels of the spots. At SMHW scales around 4° 
(10° on the sky), the data deviate from Gaussianity. The analysis is performed on all 
sky, the northern and southern hemispheres, and on four regions covering all the sky. 
A cold spot at (6 = — 57°,Z = 209°) is found to be the source of this non-Gaussian 
signature. We compare the spots of our data with 10000 Gaussian simulations, and 
conclude that only around 0.2% of them present such a cold spot. Excluding this spot, 
the remaining map is compatible with Gaussianity and even the excess of kurtosis in 
Vielva et al. 2004, is found to be due exclusively to this spot. Finally, wc study whether 
the spot causing the observed deviation from Gaussianity could be generated by sys- 
tematics or foregrounds. None of them seem to be responsible for the non-Gaussian 
detection. 

Key words: methods: data analysis - cosmic microwave background 



1 INTRODUCTION 

The Cosmic Microwave Background (CMB) is the most an- 
cient image of the universe. The measured temperature fluc- 
tuations reveal important characteristics of the primitive 
universe. Nowadays many models try to explain the primor- 
dial density fluctuations, and the key for deciding among 
them is the Gaussianity of the temperature fluctuations 
of the CMB. According to standard inflationary theories, 
these fluctuations are generated by a Gaussian, homoge- 
neous and isotropic random field, whereas non-standard in- 
flation (Linde & Mukhanov 1997, Peebles 1997, Bernardeau 
& Uzan 2002 and Acquaviva et al. 2003) and topological de- 
fect models (Turok & Spergel 1990 and Durrer 1999), predict 
non- Gaussian random fields. 

Many difficulties have to be avoided in the Gaussian- 
ity analysis, since the temperature fluctuations may come 
from different sources, which should be identified. The pri- 
mary fluctuations were originated on the last scattering 
surface (LSS), while secondary fiuctuations were produced 
in the trip that CMB photons cover from the LSS to us. 
Some examples of the latter are the fiuctuations caused by 



* e-mail: cruz@ifca.unican.es 



the reionization of the universe (Ostriker & Vishniac 1986 
and Aghanim et al. 1996), the Rees-Sciama effect (Rees & 
Sciama 1968 and Martinez-Gonzalez & Sanz 1990), or grav- 
itational lensing (Martinez-Gonzalez et al. 1997, Hu 2000 
and Goldberg & Spergel 1999). In addition we have system- 
atic effects and contaminating foregrounds, like the Galactic 
emissions or extragalactic point sources, which have to be 
considered in the Gaussianity study. 

In the last years, many groups have published anal- 
ysis of Gaussianity. The Wilkinson Microwave Anisotropy 
Probe (WMAP) 1-year aU-sky data (Bennett et al. 2003a), 
were found to be consistent with Gaussianity by the WMAP 
science team (Komatsu et al. 2003). A measure of the 
phase correlations of temperature fiuctuations (taking sev- 
eral combinations of the bispectrum into account), and the 
Minkowski functionals, were used in the latter paper. There- 
after several other groups have found non-Gaussian signa- 
tures or asymmetries in the WMAP data (Park 2003, Erik- 
sen et al. 2003, 2004a, Hansen et al. 2004a, b). Several meth- 
ods were used including Minkowski functionals, Af-point cor- 
relation functions, or local curvature methods. Also some 
other works such as Chiang et al. 2003 and Copi et al. 2003, 
detect deviations from Gaussianity, using the Internal Lin- 
ear Combination (ILC) map (Bennett et al. 2003a) and other 
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maps, derived from the WMAP data. (However, the use of 
the ILC map is not recommended by the WMAP team, be- 
cause of their complex noise and foreground properties). 

Our basic reference is Vielva et al. 2004, where a non- 
Gaussian detection in the WMAP 1-ycar data was reported. 
Convolving the data with the SMHW, an excess of kurtosis 
was found at scales axound 4°, involving a size on the sky 
close to 10°. The deviation from Gaussianity presented an 
upper tail probability of 0.4%. This deviation was localized 
in the southern hemisphere and an extremely cold spot at 
b = —57°, / = 209°, was regarded as a possible source of the 
non-Gaussianity (we call it The Spot). 

Mukhcrjec & Wang 2004 have carried out an extrcma 
analysis using the SMHW, finding an excess of cold pixels 
in all the sky at the same scales as Vielva et al. 2004. In the 
present paper, our aim is to localize the non-Gaussianity, 
specifying whether The Spot alone is the origin of the de- 
tection, and to study all the spots of the data in order to 
quantify how atypical The Spot under a Gaussian model is. 

Also recently, Larson & Wandclt 2004 have carried out 
an extrema analysis in real space, finding non-Gaussian ev- 
idences. 

The present paper is organized as follows. In §2 we de- 
scribe the data and simulations used to study the Gaussian- 
ity. The wavelet technique and the estimators are explained 
in §3. The results of our work are presented in §4. In §5 we 
discuss the possible sources of the non-Gaussian detection 
and in §6 we present the conclusions of our results. 



2 DATA AND SIMULATIONS 

We have followed a similar procedure as in Vielva et al. 
2004. Hence our data are the 1-year WMAP data, available 
in the Legacy Archive for Microwave Background Data Anal- 
ysis (LAMBDA) web page ^ , choosing for our analysis the 
combined Q-V-W map. The WMAP receivers, take their ob- 
servations at 5 frequency bands, namely K-band (22.8 GHz, 
1 receiver), Ka-Band (33.0 GHz, 1 receiver), Q-Band (40.7 
GHz, 2 receivers), V-Band (60.8 GHz, 2 receivers) and W- 
Band (93.5 GHz, 4 receivers). The combined map is a linear 
combination of the receivers where CMB is the dominant 
signal, namely Q, V, W radiometers: 



(1) 



j=3 



The Q-Band receivers are denoted with indices j = 3, 4, V- 
Band receivers with indices j = 5, 6, whereas indices j = 
7, 8, 9, 10 correspond to the four W-Band receivers. Noise 
weights Wj(x) are defined so that receivers with less noise 
have more weight: 



Wj(x) = 



Wj(x) 



Jv.,-(^) 

— ^ 



(2) 



where aoj is the noise dispersion per observation and Nj (x) 

is the number of observations performed by the receiver j 
at position x (see Bennett et al. 2003a). Foregrounds are 
corrected for each receiver, using the foreground template 

^ http://www.cmbdata.gsfc.nasa.gov/ 



described in Bennett et al. 2003b. The combined map in- 
creases the signal to noise ratio. 

The Hierarchical, Equal Area and iso-Latitude Pixeliza- 
tion (HEALPix, Gorski et al. 1999) is used in all maps. The 
resolution parameter of the original map, is Nside ~ 512 but 
we have degraded it to Nside = 256, since very small scales 
are dominated by noise. Afterwards the KpO mask was ap- 
plied in order to get rid of foreground or point source con- 
taminated pixels. This mask removes 23.2% of the pixels, es- 
pecially the strongly contaminated pixels around the Galac- 
tic plane. Finally, the residual monopole and dipole out- 
side the mask are subtracted. The mask and the foreground 
cleaned frequency maps can be found in the LAMBDA site. 

Following the same procedure as in Vielva et al. 2004, 
we have produced 10000 Gaussian simulations with the pur- 
pose of checking the Gaussianity of the data. Starting with 
the cosmological parameters estimated by the WMAP team 
(Table 1 of Spergel et al. 2003), we have calculated the Ce 
using CMBFAST (Seljak & Zaldarriaga 1996). We have gen- 
erated random Gaussian ae,n of CMB realizations and con- 
volved them with the adequate beam for each receiver. Af- 
ter transforming from harmonic to real space, uncorrelated 
Gaussian noise realizations have been added, taking into ac- 
count the number of observations per pixel {Nj{x.)) and the 
noise dispersion per observation {(Toj). At this point the sim- 
ulations were degraded to Nside = 256 and the KpO mask 
has been applied. At the end, monopole and dipole are sub- 
tracted outside the mask 

In the present work we have studied the spots above or 
below a given threshold f, therefore we have normalized the 
data and the simulations, dividing each map by its disper- 
sion, after subtracting the mean. 



3 WAVELETS AND ESTIMATORS 

Several works have been published in recent years concern- 
ing CMB studies based on wavelets. Wavelet transforms can 
be applied to separate the different components appearing 
in microwave observations (Tenorio et al. 1999, Cayon et 
al. 2000, Vielva et al. 2001a, Vielva et al. 2001b and Vielva 
et al. 2003a), in dcnoising techniques (Sanz et al. 1999a, b) 
and in Gaussianity studies (Pando et al 1998, Hobson et al. 
1999, Aghanim et al. 2003, Barreiro et al. 2000, Barreiro & 
Hobson 2001 and references listed below for the spherical 
symmetry). Wavelets increase the signal to noise ratio, al- 
lowing us to detect weak non-Gaussian signals. Furthermore 
they preserve the spatial location and the angular scale R of 
a hypothetical non-Gaussiam feature. This is a very impor- 
tant advantage since our main goal in this paper is to study 
the hot and cold spots over several thresholds. 

In our case, the ideal wavelet is the SMHW (Antoine 
J. P. & Vandcrlieynst P. 1998), which has been applied to 
non-Gaussian studies of the COBE-DMR data (Cayon et al. 
2001, 2003), to Planck simulations (Martmez-Gonzalez et al. 
2002) and to the 1-year WMAP data in our basic reference 
(Vielva et al. 2004). The SMHW optimally enhances the 
non-Gaussian signatures on the sphere (Martmez-Gonzalez 
et al. 2002). The expression of the SMHW is: 



*s(y,i?) = 



V2^N{R) 



i-Y 



,(3) 
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where R is the scale and N{R) is a normalization constant: 
N{R) = R^/l + + R^/4. The distance y on the tan- 
gent plane is related to the polar angle (6) as: y = 2tan&/2. 
Data and simulations have been convolved with the SMHW 
at different scales, to obtain the SMHW coefficient maps an- 
alyzed in this work. The study has also been performed in 
the data and simulations previous to the convolution with 
the SMHW. These maps are referred to as maps in real 
space, at scale Ro- 

We have considered the following six estimators: num- 
ber of maxima, number of minima, number of hot spots, 
number of cold spots, number of hot pixels (hot area) and 
number of cold pixels (cold area). All of them arc referred 
to a particular threshold; thus maxima, hot spots and hot 
pixels lie above a threshold v, whereas number of minima, 
cold spots and cold pixels lie below —u. The analysis has 
been done in seven regions: all sky, northern hemisphere, 
southern hemisphere, Northeast {b > 0,1 < 180), Northwest 
(6 > 0,/ > 180), Southeast {b <0,l < 180) and Southwest 
(6 < 0, / > 180). Each region has its own dispersion and has 
to be normalized separately. 



4 THE ANALYSIS 

The first step in our analysis is the convolution of data 
and simulations with the SMHW at 15 scales (i?i = 13.7, 
i?2 = 25, = 50, R4 = 75, R5 = 100, Re = 150, R7 = 200, 
Rs = 250, Rg = 300, Rio = 400, Ru = 500, R12 = 600, 
Jii3 = 750, Ri4 = 900 and R15 = 1050 arcmin). Since 
the convolution mixes up pixels of the masked region with 
unmasked pixels, the resulting map has many corrupted 
SMHW coefficients. This effect involves a loss of efficiency, 
so we applied an additional mask as explained in Vielva et 
al. 2004, extending the KpO mask to 2.5 R. The results in 
Vielva et al. 2004 were shown to be quite independent of 
the definition of this extended mask. For each map we have 
calculated the estimators on the seven previously described 
regions. The chosen thresholds are ±2.0, ±2.5, ±3.0, ±3.5, 
±4.0 and ±4.5. In the figures we have represented only ab- 
solute values for the thresholds, recalling that for cold spots 
and minima, values lie below a negative threshold, whereas 
hot spots and maxima values lie above a positive threshold. 

Once the estimators were calculated for data and sim- 
ulations, we proceeded to establish acceptance intervals at 
significance levels a (32%, 5% and 1%). Therefore we have 
sorted the 10000 values of each estimator into ascending 
numbers and excluded 10000q;/2 values from each tail of 
the distribution. The two limiting values, corresponding to 
simulations 1 + lOOOOa/2 and 10000(1 - a/2) define the ac- 
ceptance interval containing a probability of 1 — a. The re- 
maining probability is a/2 above and below the interval. In 
all figures this acceptance intervals are plotted, the 32% in- 
terval corresponds to the inner band, the 5% interval to the 
middle band and the 1% significance level, to the outer one. 

The results for the different regions are presented in the 
next three subsections. 

4.1 All sky, North and South 

Deviations from Gaussianity have been detected at scales 
i?B, Ra and Rg. For threshold 3.0, the results are repre- 



sented in Figure 1. At scale R5 the hot area lies outside the 
99% acceptance interval. However we have studied this case 
carefully, noting that this threshold is the only one were 
data lie outside the intervals, whereeis the number of hot 
spots and number of maxima do not depart from Gaussian- 
ity, even at this scale. Dividing the sky into two hemispheres 
we checked that this was not a localized effect either. Even 
contiguous scales were compatible with Gaussianity, hence 
we concluded that this non-Gaussian feature was not signif- 
icant. 

On the contrary, in the case of cold areas, deviations are 
observed at scales Ra and Rg at several thresholds as can be 
seen in Figure 2. We represent in Figure 2 the number of 
cold spots, cold area and number of minima, for all consid- 
ered thresholds at scales Rg and Rg, in order to observe what 
is happening there with more detail. The number of minima 
and cold spots was exactly one at thresholds 4.0 and 4.5, 
reaching the borderline of our 99% acceptance interval. For 
the number of cold spots, at threshold 4.0 for scale Rg and 
thresholds 3.5 and 4.0 for scale Rg, the 95% and 99% accep- 
tance intervals coincide. This can happen when the number 
of spots is very low. Consider for example a binomial dis- 
tribution where the rmmbcr of spots can only take values 
or 1. The way we define the acceptance intervals determines 
that for such cases some of these intervals must coincide. 
The most striking non-Gaussian signature was found for the 
cold area at scales -Rs and -R9 where the data exhibited an 
extremely high immber of cold pixels at thresholds over 3. 
The cold area lies outside the 99% acceptance interval at 
thresholds above 3.0 in the two scales presented in Figure 
2. 

At thresholds 4.0 and 4.5 the mentioned observations 
reveal that the non-Gaussianity is only due to one par- 
ticular spot, which reaches a minimum value of — 4.7cr at 
(6 = —57°, I = 209°) and scale Rg. At lower thresholds sev- 
eral spots contribute to the observed deviation. A precise 
analysis of The Spot is presented in the following sections. 
The data suggest that we are dealing with a very cold and 
big spot. These results agree with the results reported in 
Vielva et al. 2004, since the non-Gaussianity has been de- 
tected at the same scales. The hot spots did not show any 
non-Gaussian evidence. Furthermore, subtracting the cold 
pixels from the hot pixels, a strong hot-cold asymmetry was 
revealed at scales Ra, Rg and Rio, see Figure 3. The lower 
tail probabilities displayed in Table 1 show non-Gaussian 
values of up to 99.85%. 

Our purpose was to locate the non-Gaussian sources. 
We studied the northern and southern hemispheres sepa- 
rately, expecting to find non-Gaussian results in the south- 
ern hemisphere because The Spot is located there. Results 
presented in Figure 4 show that the northern hemisphere 
is compatible with the Gaussian simulations whereas the 
southern hemisphere shows a clear deviation from the ac- 
ceptance intervals. To show the North-South asymmetry, the 
number of hot and cold pixels in the South were subtracted 
from the Northern ones. Results are presented in Figure 5. 
Deviations from Gaussianity can be clearly observed in the 
number of cold pixels at thresholds 3.0 and above. Previ- 
ous works such as Vielva et al. 2004, Erikscn et al. 2004a or 
Hansen 2004b have already reported North-South asymme- 
tries, which axe confirmed here. 

The marked observed asymmetry, prompted us to di- 
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All sky, threshold = 3.0 



All sky, threshold = 3.0 




All sky, threshold = 3.0 




All sky, threshold = 3.0 




15 



All sky, threshold = 3.0 



All sky, threshold = 3.0 





Figure 1. This figure shows from left to right and top to bottom, number of maxima, number of minima, number of hot spots, number 
of cold spots, hot and cold area for each of the considered scales at threshold 3.0. Scale Rq corresponds to real space. The combined 
map values are plotted as asterisks. The acceptance intervals for the 32% (inner), 5% (middle) and 1% (outer) significance levels, given 
by the 10000 simulations are also plotted. Non-Gaussianity is found at scales in the hot area, Rg and Rg in the cold area. At other 
scales the data always lie at least within the 99% acceptance band. For scale R^ the deviation occurs only for this threshold and scale 
whereas scales Rg and Rg show a more significant detection. 
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Scale = R , all sky 



Scale = R all sky 





4.5 2 



thresholds 
Scale = R , all sky 



Scale = R , all sky 




thresholds 
Scale = R , all sky 



Scale = R , all sky 




Figure 2. Number of minima, number of cold spots and cold area are shown for scale _Rs at the left column and scale Rg at the right 
column. As in figure 1 the asterisks represent the combined WMAP data, and the acceptance intervals for the 32% (inner), 5% (middle) 
and 1% (outer) significance levels, given by the 10000 simulations are also shown. The data lie outside the acceptance intervals for the 
area, at thresholds over 3.0. The borderline of the 5% and 1% acceptance intervals coincide at some high thresholds for the number of 
spots, because at this scales the number of spots is very low. 
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Scale = R , hot-cold Scale = R , hot-cold 




2 2.5 3 3.5 4 4.5 2 2.5 3 3.5 4 4.5 

thresholds thresholds 
Scale = R^^ , hot-cold 



8000 
6000 
4000 
2000 
o 

I 

-2000 
-4000 
-6000 



2 2.5 3 3.5 4 4.5 

thresholds 

Figure 3. Hot-cold asymmetry, in the area. The excess of cold pixels at scales ijg, Rg and Rio happens at high thresholds. Data and 
acceptance intervals are represented in the same manner as in previous figures. 




Scale 


threshold 


probability 


Rs 


3.0 


0.37% 




3.5 


0.26% 




4.0 


0.44% 




4.5 


0.65% 


R9 


3.0 


0.39% 




3.5 


0.15% 




4.0 


0.19% 




4.5 


0.22% 


Rio 


3.0 


18.22% 




3.5 


1.04% 




4.0 


0.48% 



vide the sky into four regions with the purpose of studying 
a Southeast-Southwest asymmetry. 



Table 1. Lower tail probabilities of having the hot-cold asym- 
metry (in the number of pixels) of our data under a Gaussian 
model, at different scales and thresholds. Most of the values are 
below 1%. At scale Rio and threshold 4.5 the number of hot and 
cold pixels are both zero, because the threshold is too high and 
therefore we do not show this threshold for the mentioned scale. 



4.2 Four regions 

We analyzed four regions independently by splitting each 
hemisphere into two parts, I < 180 (East) and I > 180 
(West). The results presented in Figure 6 clearly reveal the 
location of the non-Gaussian signatures. The Southwest is 
the only region where the data lie outside the acceptance in- 
tervals. As expected, the region containing The Spot is not 
compatible with Gaussianity, whereas the other three re- 
gions are compatible with a Gaussian behaviour. Note that 
although Eriksen et al. 2004b found the ILC weights to have 
particular values in this region, this does not affect our re- 
sults since we use the combined, cleaned Q-V-W map. 

In the next section we quantify the significance of The 
Spot. 
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Scale = R„ , north 




3 3.5 

thresholds 
Scale = R„ , nortti 




3 3.5 
thresholds 



Scale = Rg , soutti 




3 3.5 

thresholds 



Figure 4. Cold area in the northern and southern hemispheres at scales and Rg. The southern hemisphere shows a non-Gaussian 
behaviour. Again we represent data and acceptance intervals as in Figure 1. 



4.3 The Spot 

Our aim in this section is to quantify the probability of find- 
ing a spot like the one found at {b = —57°, I — 209°) in a 
Gaussian, homogeneous and isotropic random field, and to 
check whether the map without The Spot is compatible with 
Gaussianity or not. Comparing The Spot with the biggest 
cold spot of each simulation we can estimate the probability 
of observing such a spot for the Gaussian model. We show 
the histograms with the biggest cold spot of each simulation 
at threshold 3.0 and scales Rs and Rg in Figure 7. The 
results for all thresholds are summarized in Table 2. Note 
that some simulations do not have any spots at high thresh- 
olds and hence they do not appear in the histograms but 



we take them into account to estimate the probabilities. All 
probabilities are below 0.7%. The lowest value is 0.18% and 
implies a non-Gaussian detection at the 99.82% level. 

At this point we can make the hypothesis that the data 
could be explained as the sum of a Gaussian, homogeneous 
and isotropic random field, plus a non-Gaussian spot which 
is not generated by this field. With the purpose of checking 
our hypothesis, we have compared the cold area of data and 
simulations, at scales Rs, Rg and thresholds 3.0, 3.5 where 
the data present more than one spot. First we have esti- 
mated the upper tail probabilities of finding the total cold 
area (including all cold spots), counting how many simula- 
tions present a greater or equal cold area, obtaining very 
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Scale = R„ , north-south 




3 3.5 
thresholds 




3 3 5 

thresholds 



Figure 5. The two upper panels represent the North-South asymmetry for hot and cold area, at scale Rg and the two pannels below, 
the mentioned asymmetries at scale Rg. The southern hemisphere shows an excess of cold pixels. Acceptance intervals and data are 
represented as in previous figures. 



low values (see Table 3). If we subtract The Spot from the 
cold area in the data, and calculate again the probabilities, 
we appreciate how the remaining area is compatible with 
Gaussianity, since the upper tail probabilities grow in a fac- 
tor of 100. This result is represented in Figure 8. The line 
on the right hand side is the total cold area of the data 
whereas the one on the left hand side is the area remaining 
after subtracting The Spot. The increase in probability can 
easily be appreciated. Here we should have considered that 
we were using fewer pixels in the data than in the simu- 
lations because The Spot has been subtracted only in the 
data. But since these pixels represent only about 0.5% of 
the total pixels at scales -Rs and Rg, we can neglect them, 
without modifying substantially the results. Hence regard- 



ing The Spot as a non-Gaussian outlier, the remaining data 
are compatible with Gaussianity. 

To finish this section, we want to remark some charac- 
teristics of The Spot. Before convolving the combined map. 
The Spot appears as several smaller, resolved spots with a 
minimum temperature of —398fiK at {b — —56°,/ — 210°). 
The average value of the four most prominent spots, listed 
in table 4 is ~34:8^K. However, these four spots are not 
exclusively responsible for The Spot, since removing them, 
the mean value of the remaining pixels forming The Spot, 
is close to —la. We have filtered the combined map with a 
Gaussian of scale 4°, the characteristic scale given by the 
SMHW analysis. The Spot appears with an amplitude of 
-78nK at (b = -56°, Z = 210°) as summarized in Table 4. 
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Scale = R„ , northwest 




Scale = R , northeast 



3 3.5 
thresholds 
Scale = R„ , southwest 





Scale = Rg , northwest 



Scale = R„ , northeast 




2 2.5 3 3.5 4 4.5 
thresholds 
Scale = R„ , southwest 



2.5 3 3.5 4 

thresholds 
Scale = R„ , southeast 




3 3.5 4 
thresholds 



Figure 6. Cold area in four regions, scales i?g and ijg. The only non-Gaussian region is the southwest. The acceptance intervals are 
plotted as in Figure 1. 



4.4 Skewness and kurtosis 

At scales Re, and Rq, the kurtosis presented non-Gaussian 
values in Vielva et al. 2004. We have repeated the analy- 
sis of Vielva et al. 2004, masking The Spot in the wavelet 
coefficient maps, to show to which extent The Spot is re- 
sponsible for the excess of kurtosis. The masked pixels, were 
the pixels of The Spot lying above threshold 3.0 at scale 
Rs. The results are presented in Figure 9. The stars repre- 
sent the data and the circles, the data without the Spot. By 
subtracting The Spot, the decrease of the kurtosis is clearly 
observed, being now compatible with Gaussianity. The ac- 
ceptance intervals are the same as in Vielva et al. 2004, and 



they were not recalculated masking the pixels of The Spot in 
the simulations, since the number of masked pixels is negli- 
gible with respect to the total number of pixels. However the 
decrease in the kurtosis is so huge, that a slight modification 
of the intervals would not affect our conclusions. The skew- 
ness is still compatible with Gaussianity. Note that the Spot 
was masked after convolving with the wavelets. Masking The 
Spot before convolving, the decrement of the kurtosis is even 
higher. This results show again the Gaussian behaviour of 
the data without The Spot. Hence we can conclude that the 
excess of kurtosis is exclusively due to The Spot. 
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Scale = R , thresliold=3.0 , all Sky 



Scale = R , threshold=3.0 , all sky 




500 1000 
number of pixels 



1000 1200 1400 1600 1800 2000 
number ot pixels 



Figure 7. Histogram of all biggest cold spots, threshold 3.0, scales and Rg. The vertical line represents The Spot. 



Scale = R , threshold = 3.0 , all Sky 



Scale = R , threshold = 3.0 , all sky 




Figure 8. Histogram of all cold areas, threshold 3.0, scales Rs and Rg. The data without The Spot (line on the left) are compatible 
with Gaussianity, whereas the data with The Spot (line on the right), are non-Gaussian with probabilities 99.82% and 99.61% at this 
threshold. 



Scale 


threshold 


probability 


R& 


3.0 


0.34% 




3.5 


0.32% 




4.0 


0.41% 




4.5 


0.65% 


R9 


3.0 


0.38% 




3.5 


0.21% 




4.0 


0.18% 




4.5 


0.22% 



Table 2. Upper tail probabilities of having The Spot under a 
Gaussian model, at different scales and thresholds. All probabili- 
ties are smaller than 1%. 



5 SOURCES OF NON-GAUSSIANITY 

Although Vielva et al. 2004 showed that systematics, fore- 
grounds and variations of the power spectrum were not re- 
sponsible for the non-Gaussian effect shown in the kurtosis, 
we wanted to check again their influence in the non-Gaussian 
resuhs obtained in the present analysis. 



5.1 Systematics 

First we studied the effects of systematics related to instru- 
mental features (noise and beam), generating four sets of 10 
simulations. The first two sets are normal simulations with 
noise and beams. In the third set we made the same simula- 
tions as in the first set but without noise, and in the fourth 
set we took the simulations of the second but without beam. 



© 0000 RAS, MNRAS 000, 000-000 



A non-Gaussian Spot in WMAP 11 



Skewness Kurtosis 




5 10 15 5 10 15 

#R #R 



Figure 9. This figure shows the values of skewness (left) and kurtosis (right) for all the considered scales. The combined map values 
are plotted as stars and as written in Vielva et al. 2004, at scales ijg and ijg they lie outside the three acceptance intervals, which are 
represented as in previous figures. Masking those pixels of The Spot which are above threshold 3.0 at scale R^, we obtain the results 
represented by circles. The huge decrement of the kurtosis makes the data compatible with Gaussianity. 



Scale 
Hg 


threshold 
3.0 
3.5 
4.0 
4.5 


P with Spot 
0.18% 
0.28% 
0.45% 
0.65% 


P without Spot 
14.79% 
18.28% 




3.0 


0.39% 


30.53% 




3.5 


0.18% 


17.68% 




4.0 


0.19% 






4.5 


0.22% 





Table 3. Upper tail probabilities of having the cold area mea- 
sured in the data, under a Gaussian model The third column 
displays the probabilities considering The Spot and the right col- 
umn shows the probabilities subtracting The Spot from the data. 
At thresholds 4.0 an 4.5 only The Spot is present, hence it makes 
no sense considering probabilities without The Spot. 



Combined WMAP 



Amplitude (p,K) Position (b, I) area {y = —2) 

-346 (-60°, 213°) 46 

-398 (-56°, 210°) 67 

-331 (-54°, 211°) 42 

-317 (-56°, 203°) 88 



Gaussian (4°) 

Amplitude (fiK) Position (b, I) area {y = —2) 

-73 (-57°, 209°) 699 



Table 4. Before filtering the combined map. The Spot appears 
resolved in several spots. At the top we show three characteris- 
tic values of the most prominent spots, namely the area below 
threshold -2, amplitude and position of the minimum. Filtering 
with a Gaussian of 4° all these spots are convolved, becoming 
one big spot. The characteristics of this spot are shown at the 
bottom. 



Comparing the first and the third set, we can see how the 
noise afltects the number of pixels in a spot or the total cold 
area, and comparing sets two and four, we check the influ- 
ence of the beams. We have compared the spots and the 
total area of these sets at scales J?g and i^g, and threshold 
3.0. In all cases the mean relative variation of the area of 
one spot was around 2% and the mean relative variation of 
the total area in a simulation around 1%. This values are 
negligible and confirm that noise and beams do not play a 
signiflcant role in our results, although the number of con- 
sidered simulations is not very high. The previous results of 
Vielva et al. 2004 support our conclusions. 

Another possible source of non-Gaussianity could be the 
influence of any rare receiver. We have analyzed the spots 
detected by the 8 Q-V-W receivers, Ql, Q2, VI, V2, Wl, 
W2, W3 and W4, independently. The results for scales JJg 
and 7?9 are plotted in Figure 10. Although W2 detects less 
pixels than the other receivers, all of them detect The Spot, 
and are close to the line representing the values of the com- 
bined map. Hence we verify that our detection is not due to 
any deficient receiver. 

5.2 Foregrounds 

Non-Gaussianity can be generated by foregrounds due to 
synchrotron, free-free and thermal dust emissions. All fore- 
grounds show a clear frequency dependence and hence if 
our Spot is generated by foregrounds its area should also 
be frequency dependent. Therefore, in Figure 11 we have 
compared the area of The Spot for each channel, namely 
K(22.8 GHz), Ka(33.0 GHz), Q(40.7 GHz), V(60.8 GHz) 
and W(93.5 GHz) at scales Rg, and -Rg and thresholds over 
3. The horizontal lines denote the combined map values, 
whereas the other symbols correspond to the area of The 
Spot at different thresholds for each channel. 

The first two channels are not foreground corrected and 
so they may not match the results of the other channels 
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Figure 10. Area of The Spot for all receivers, at threshold 3.0 and scales Rs and Rg. The line represents the combined map values 
whereas the other symbols denote the diflferent receivers. The Spot is detected in all receivers, hence no rare receiver is generating the 
non-Gaussian signal. 



400-, * 





threshold 3.0 





threshold 3.5 




threshold 4.0 


□ 


threshold 4.5 


o 


K-2.65Ka 




combined map 




* 


threshold 3.0 





threshold 3.5 




threshold 4.0 


□ 


threshold 4.5 


o 


K-2.65Ka 




combined map 



Figure 11. Frequency dependence of the area of The Spot, at thresholds over 3.0 and scales Rg and Rq. The lines represent the combined 
map values. From top to bottom thresholds 3.0 (asterisks), 3.5 (diamonds), 4.0 (stars) and 4.5 (squares) are plotted. The circles are the 
values for the map K — 2.&bKa, where the synchrotron emission should be canceled. These circles are close to the combined map values, 
hence synchrotron emission is not the cause of The Spot. 



and of the combined map. In fact £is can be seen in Figure 
11 both channels deviate slightly from the horizontal lines, 
representing the combined map values. Since these channels 
are not foreground corrected, we can attribute the deviation 
to the synchrotron radiation, which dominates at these fre- 
quencies. The synchrotron emission is expected to grow a 
factor 2.65 from 33 GHz to 23 GHz,^ hence if we subtract 
2.65 times the Kar-map from the K one, we get rid of the 

■t A power law is assumed for the frequency dependence of the 
synchrotron emission: Tsyn{i^) <x 7sj;„(i/o) I i^/i'o ) i as pro- 
posed by Bennett et al. (2003b). 



contaminating emission. Considering Figure 11, we confirm 
that the circles corresponding to K-2.65Ka are very close to 
the combined map values. 

These results support the conclusions reached in Vielva 
et al. 2004 regarding the influence of foregrounds in the non- 
Gaussian detection The independence of the amplitude of 
The Spot with frequency, was already shown in the men- 
tioned paper. 

5.3 Power spectrum dependence 

Since several anomalies and asymmetries have been found 
in the low multipoles of the power spectrum, we should dis- 
cuss to which extent the result depends on the power spec- 
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trum. We have used the best fit WMAP power spectrum to 
perform the 10000 Gaussian simulations, but the uncertain- 
ties in the cosmological parameters and hence in the power 
spectrum could affect our results. We have performed three 
sets of 50 simulations. The first set corresponds to the best 
fit WMAP power spectrum. The other two sets were gen- 
erated with power spectra differing by ±lcr from the best 
fit WMAP power spectrum. One corresponds to the 'lower 
limit' power spectrum, obtained subtracting the Icr error 
estimated by the WMAP team, and the other one to the 
'upper limit' power spectrum, adding the mentioned la er- 
ror to the best fit spectrum. Comparing the first set with 
the two others, we can study how the spots are affected by 
the choice of different power spectra. 

At scales Rs and R\), and threshold 3.0 the mean rela- 
tive variation of the area of a particular cold spot is much 
lower than 1%. Therefore, after these negligible variations, 
the significance of The Spot remains unchanged. Also Vielva 
et al. 2004 found a negligible power spectrum dependence of 
the acceptance intervals for the kurtosis. Hence the choice 
of different power spectra does not affect significantly our 
results. 



5.4 Intrinsic anisotropies 

Once we have discarded systematics and foregrounds as the 
cause of our detection, other sources have to be consid- 
ered. For instance the Sunyaev-Zeldovich effect could pro- 
duce a cold spot. This effect occurs in clusters, when high 
energetic electrons collide with CMB photons originating a 
decrement of the temperature in our range of frequencies. 
Here arise two problems, namely the angular scale and the 
amplitude of The Spot. The non-Gaussianity is found at 
scales around 4° implying a size of about 10° on the sky. 
Observing in real space this region. The Spot appears re- 
solved in several smaller very cold spots with a minimum 
temperature around — 398/iK. We have looked for any ex- 
tragalactic object which could cover this angular scale at 
coordinates near (6 = — 57°,Z = 209°). We found a group of 
galaxies belonging to the local supercluster at a distance of 
about 20 Mpc subtending a similar angle on the sky. Some 
of these galaxies match the resolved spots observed in the 
data. However the mass and temperature of the gas nec- 
essary to reach the amplitude of our Spot, are similar to 
the values found in rich clusters and therefore much higher 
than the amounts estimated for groups of galaxies (see Tay- 
lor et al. 2003). We would need a big cluster, such as the 
Coma one, to reach such an amplitude and it should be near 
enough to cover 10 degrees on the sky. In the neighborhood 
of (6 = —57°, I = 209°) no such object is found. This is in 
agreement with the WMAP results on foregrounds (Bennett 
et al 2003b) where the Sunyaev-Zeldovich effect was found 
to be negligible except for the most prominent nearby clus- 
ter, Coma, observed with a signal to noise ratio of ~ 2. Even 
more, we have also considered the ACO catalogue and the 
All-Sky ROSAT maps (Snowden et al. 1997) at 0.1, 1.2 and 
2 keV and neither any ACO cluster nor any special X-ray 
emission was found at the The Spot position. However, the 
ROSAT maps present some particular problems: the bright- 
est point sources, as well as a large fraction of clusters, have 
been removed from it. In addition, some small fractions of 



the sky were not observed and, unfortunately, one of them 
is very close to our object. 

Another possible source is the Rees-Sciama effect, 
(Martinez-Gonzalez et al. 1990, Martinez-Gonzalez & Sanz 
1990). An extremely massive and distant superstructure 
would be a clear candidate to cause a cold and big secondary 
anisotropy. Even topological defects like global monopoles or 
textures (Turok & Spergel 1990) could have cooled the CMB 
photons, to produce such spot. Cosmic strings have charac- 
teristic scales around arcminutes, hence we do not expect 
them to be behind this non-Gaussian detection. 



6 CONCLUSIONS 

Motivated by the non-Gaussianity found in the WMAP 1- 
year data using the SMHW, we have performed an analysis 
of the spots in the SMHW coefficients map, aimed to locate 
possible contributors in the sky. An extremely cold and big 
spot is detected. This spot ( The Spot) is seen in the SMHW 
coefficients at scales around 4° (implying a size of around 
10° on the sky) and at Galactic coordinates b = —57°, I = 
209°. The probability of having such spot for a Gaussian 
model is of only w 0.2%, which implies that, if intrinsic. 
The Spot has not been originated by primary anisotropies in 
the standard scenario of structure formation since standard 
inflation predicts Gaussian fluctuations in the matter energy 
density and therefore in the CMB temperature fluctuations. 
When this spot is not considered in the analysis the rest of 
the data seem to be consistent with Gaussianity. 

In order to identify the source of The Spot we have per- 
formed several tests related to systematic effects and fore- 
grounds. We have checked that uncertainties in the noise or 
in the beam response have a negligible effect in our results at 
the relevant wavelet scales. Looking at the maps correspond- 
ing to the different receivers, we sec a clear consistency in 
the area, amplitude and position of The Spot. Hence our de- 
tection is not due to any deficient receiver. In relation to the 
possible foregrounds contribution, we have looked for possi- 
ble frequency dependences in the amplitude and area of The 
Spot. Again both quantities show a nice consistency with a 
constant line in the range from 23 to 94 GHz. Whereas the 
Galactic foregrounds show a very different frequency depen- 
dence with respect to the constant behaviour, the SZ effect 
does not separate much from it in that frequency range. 
A comparable spot could be produced either by the Coma 
cluster at a much closer distance, or by several rich clusters 
at the actual distance of Coma. We have checked that no 
nearby rich cluster of galaxies is located in the position of 
The Spot. 

Finally, intrinsic fluctuations cannot be rejected as the 
source of The Spot. In particular, a massive and distant 
super-structure could in principle produce a decrement as 
the one observed through the Rees-Sciama effect (Martinez- 
Gonzalez and Sanz 1990). This massive structure (of order 
of at least 10^® Mq) should be placed far enough because 
otherwise it would have been detected previously. Alterna- 
tively, more speculative possibilities arc tojjological defects 
(monopoles or textures) or non-standard inflationary sce- 
narios. Even more, a combination of secondary anisotropies 
with primary ones, cannot be rejected as the source of our 
non-Gaussian spot. For instance a possibility could be a 
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combination of the Sunyaev-Zeldovich effect with a Sachs- 
Wolfe plateau. 
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